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RANDOM  COVERAGE  OF  A  CIRCLE  WITH  APPLICATION 


TO  A  SHADOWING  PROBLEM 


1.  INTRODUCTION 

The  problem  under  consideration  in  the  present  study  belongs  to  the  class  of 
problems  of  random  coverage  of  a  circle  by  randomly  placed  arcs  having  random 
length.  The  specific  problem  studied  here  is  motivated  by  a  shadowing  problem, 
according  to  which  the  random  arcs  on  the  circle  are  shadows  cast  by  randomly 
scattered  disks  on  the  plane.  The  model  assumes  that  the  centers  of  the  disks 
form  a  homogeneous  Poisson  process  on  the  plan,  with  intensity  v  per  unit  area. 
Furthermore,  given  the  number  of  disks  centered  in  a  specific  Borel  set  in  the 
plane,  the  model  assumes  that  their  diameters  are  identically  distributed  inde¬ 
pendent  (i.i.d)  random  variables  with  a  known  distribution.  Assuming  that  the 
circle  is  not  intersected  by  any  disk  and  its  center  (the  source  of  light)  is 
uncovered,  shadow  arcs  on  the  circle  are  defined  as  the  central  projections  of 
random  disks  which  lie  entirely  within  the  set  inscribed  by  the  circle.  Accord¬ 
ingly,  the  model  assumes  that  the  number  of  shadow  arcs  covering  the  circle  is  a 
Poisson  random  variable.  Given  this  number,  the  centers  of  the  disks  are  condi¬ 
tionally  independent  having  a  uniform  distribution.  The  length  of  the  shadow 
arcs  are  conditionally  i.i.d.  random  variables,  having  a  common  distribution  on 
[o,  it].  The  main  objective  of  the  present  study  is  to  obtain  the  distribution 
of  a  measure  of  vacancy  of  arcs  on  the  circle.  The  measure  of  vacancy  of  an  arc 
is  the  total  length  of  the  portion  of  the  arc  which  is  in  the  light.  Equations 
for  the  moments  and  the  moment  generating  function  of  the  measure  of  vacancy  of 
arcs  of  length  t,  0  <  t  <  2*,  and  given  in  section  4.  Furthermore,  for  0  < 
t  <  n  a  formula  of  the  Laplace  transform  of  this  moment  generating  function  is 
developed.  The  derivations  in  Section  4  are  based  on  formulae  of  vacancy  proba¬ 
bilities  derived  in  Sections  2  and  3.  More  specifically,  in  Section  2  we  pro¬ 
vide  general  formulae  for:  (1)  the  probability  that  an  arc  of  length  t  on  the 
circle  is  not  covered;  (2)  the  probability  that  a  finite  set  of  specified  points 
are  simultaneously  uncovered.  These  formulae  are  further  developed  in  Section  3 
in  terms  of  the  stochastic  specifications  of  the  random  arcs.  In  Section  5  the 
results  are  applied  to  the  particular  shadowing  problem  under  consideration. 
Analytic  formulae  for  the  numerical  determination  of  the  moments  of  any  order, 
for  a  particular  example,  are  given  in  Section  6. 

The  literature  on  the  coverage  problem  is  very  extensive.  Robbins  [4]  de¬ 
rived  the  moments  of  the  total  coverage  of  an  interval  on  the  line  by  random 
segments  of  fixed  size.  These  results  were  later  extended  and  generalized  by 
Robbins  [5],  Domb  [2],  and  Takacs  [9].  Other  related  results  are  presented  by 
Solomon  [8j.  Siegel  presented  in  [6]  moments  of  the  measure  of  vacancy  of  the 
circle  when  the  coverage  is  by  a  fixed  number  of  random  arcs  having  random 
length.  In  a  following  paper  [7]  Siegel  provided  formulae  for  the  moments  and 
the  distribution  of  the  measure  vacancy  on  a  circle,  which  is  covered  by  a  given 
number  of  random  arcs  of  fixed  length.  The  motivating  shadowing  problem  led  us 
to  further  developments  over  those  of  the  previous  papers,  although  the  basic 
approach  to  the  evaluation  of  moments  is  essentially  the  same  as  that  of 
Robbins. 
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The  shadowing  problem  did  not  receive  much  attention  in  the  literature. 
Chernoff  and  Daly  [1]  considered  a  similar  shadowing  problem  when  the  shadows  of 
random  disks  are  cast  on  a  straight  line.  They  provide  the  methodology  for  de¬ 
veloping  the  distributions  of  the  length  of  intervals  which  are  entirely  in  the 
light  or  entirely  in  shadow.  It  should  be  remarked  that  the  shadowing  process 
discussed  in  the  present  paper  is  identical  over  the  interval  [o,  with  an 
M/G/«  queuing  process. 

2.  THE  COVERAGE  MODEL  AND  FUNDAMENTAL  RESULTS 

Consider  a  circle,  C,  of  radius  one  centered  at  the  origin.  Let  A},  A?,..., 
An  be  N  arcs  placed  at  random  on  C.  N  is  a  random  variable  having  a  Poisson 
distribution  with  mean  2uX.  Given  that  N  =  n,  the  centers  of  Aj,...  An  are  con¬ 
ditionally  independent  and  uniformly  distributed  on  C.  Let  Xi  (i  *  l,...n)  de¬ 
note  the  arc  length  of  A-j.  It  is  assumed  that  Xi...,  Xn  are  i.i.d.  random  vari¬ 
ables  having  a  c.d.f.  F(x)  on  [o,  *]. 

It  is  well  known  that  the  probability  of  covering  any  specific  point  on  C  by 
a  randomly  placed  arc  is  y  *=  E{X}/2*.  This  result  will  be  obtained  as  a  special 
case  of  a  more  general  result  derived  in  the  present  paper.  Let  PT  designate  a 
point  on  C  having  polar  coordinates  (1,  t),  0  <t<2u.  Let  Q(t)  denote  the  number 
of  arcs  which  cover  PT.  Given  N  e  n,  the  conditional  distribution  of  Q(t)  is 
the  binomial  B(n,Y).  Accordingly  the  (total)  distribution  of  Q(t)  is  the  Pois¬ 
son  with  mean  p  ■  2nVr  *  XE{X).  Notice  that  due  to  the  symmetry  in  the  model 
the  distribution  of  Q(t)  does  not  depend  on  t. 

Consider  a  specified  arc  on  C  of  length  t,  o<t<2»,  connecting  the  points  Ps 
and  Pe+t*  Let  qn(t)  denote  the  conditional  probability  that  such  an  arc  is  com¬ 
pletely  uncovered  by  random  arcs,  given  N  *  n.  This  conditional  probability  is 


q„(t)  -  P{  V  Q ( t )  =  0|N*n} ,  (2.1) 

S<T<S+t 


where  VQ(t)  denotes  the  maximum  of  Q(t)  over  the  specified  interval.  Since 

T 

the  random  arcs  are  conditionally  independent,  given  N=n,  qn(t)  *  (qj  (t))n,  for 

all  t,  o<t<2*.  Accordingly,  the  probability  that  a  specified  arc  on  C  of  length 
is  uncovered  is 


Explicit  derivation  of  q^t)  is  given  in  Section  3.  In  particular,  for  t«0,  we 

obtain  that  qx (O)-l--r  and  the  probability  that  ar\y  given  point  on  C  is  uncovered 
is 


£  *  q(0)  ■  e"p  (2.3) 


Let  Ps  ,...,  Pg  be  r  specified  points  on  C,  where  r>2  and  0<Sj  <s2<...<sr<2». 
r 

Let  tfta  2*-sr+s,  ,  ti»s  -s  ,  i*l . r-1,  be  the  length  of  the  arcs  connecting 

0  r  1  i+1  i 

these  points.  The  conditional  probability,  given  N*n,  that  all  these  r  points 
are  uncovered  is  function  of  (t( .  tr-l)  defined  by 

pn  (t,,...,  tr.i)  -  P  {  V  Q(s-j )  *  OjN*n}  (2.4) 

i*l,...r 


r-1 

0<ti  r-1)  and  l  ti<2*. 

i-1 


Again,  due  to  the  conditional  independence,  pn(tj,...,  tr.j)  *  (Pjftj,...,  tr.j))n 

r-1 

for  all  (t,,...,  tr_i).  Notice  that  t0  «  2*  -  l  t,.  Finally,  the  (total)  prob- 
1  i«l 

ability  that  the  r  points  are  uncovered  is 

p(tj . tr_i)  «  exp  {-2*x(l-p1(t1,...,tr.i))}.  (2.5) 

Explicit  formula  for  p^tj,...,  tr_i)  will  be  given  in  Section  3. 
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3.  VACANCY  PROBABILITIES 

In  the  present  section  \ne  develop  explicit  formulae  for  qj(t)  and  Pj(tj,..., 
tr_i).  As  defined  in  the  previous  section,  qj(t)  is  the  probability  that  a  ran¬ 
domly  placed  arc  does  not  intersect  a  specified  arc  of  length  t.  This  probability 
is  given  by 

i  2»-t 

q,(t)  •  —  /  F(x)dx.  (3.1) 

1  2*  o 

Introduce  the  auxiliary  function 

2* 

♦(t)  »  /  [1  -  F(x)]dx,  0<t<2».  (3.2) 

t 

Notice  that  4>(0)  «  E (X)  and  that,  ♦(t)«0  for  all  t>*.  Accordingly, 

qx(t)  -  1  -  L  [t  4  E (X)  -  ♦(2*-t)],  0<t<2*.  (3.3) 

From  (3.3)  one  obtains  the  previously  mentioned  result  that  9,(0)  *  l-E{X}/2*  *  1-y 
and 

q(t)  »  ?exp  {-X(t-4>(2*-t))} ,  0<t<2*.  (3.4) 
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Another  important  result  that  can  be  obtained  from  (3.4)  is  the  distribution  of 
the  length  of  an  uncovered  arc  starting  at  an  uncovered  point  Ps.  More  specifi¬ 
cally,  consider  the  r.v. 


Hs  ■  sup  (t;  V  Q(t)  ■  0,  0<t<2*} 

S<T<S+t 


(3.5) 


where  Q(s)  *  0.  It  follows  that 


P  (Hs  >  t)  -  q(t)/C 
,  e’xt  ,0<t<* 

*  g  -Xt+X$(2tr-t )  ,ir<t<2* 

0  ,2w<t 

Notice  that  the  distribution  of  Hs  has  a  jump  point  at  t«2»  and 

P  (He  -  2*}  «  e-2llX/Z. 


(3.6) 


(3.7) 


For  the  derivation  of  an  explicit  formula  for  px (tj ,. . . ,  tr.i),  let  B0  denote  the 

event  that  a  single  random  arc  lies  entirely  within  the  arc  between  Ps  and  Ps 

r  l 

and  Bi  (i«l .  r-1)  the  event  that  the  random  arc  lies  entirely  between  Ps.  and 


Ps  .  Thus, 
i+1 


P  {Bj}  -  q(2w  -  ti) 


‘-t-  [tj  +  4»(ti )  -  E{X}],  1*0,.*.,  r-1# 


(3.8) 


Finally,  since  Br_i  are  mutually  exclusive  and  their  union  is  the  event 

{.  .  V  Q(Si)  »  0}, 


r-1 


tr«i)  *  I  PfB-j) 
1  i»o 


r  i  r-1  ,  . 

-  1  -  -  E(X)  ♦  —  I  *(ti). 

2»  2»  i »o 


(3.9) 


Substitution  of  (3.9)  in  (2.5)  yields 


r-1 

p(t.,...,  tP_i)  -  Cr  exp  {X  l  ♦(ti)}. 
1  i«o 


(3.10) 


In  particular,  for  r«2,  the  joint  probability  that  two  specified  points,  Ps  and 
Ps+t,  are  uncovered  by  random  arcs  is 


p(t) 


£2eX*(t)  ,0<t<» 

C2  .t-» 

52ex^(2*-t)  ,T<t<2ir  . 


(3.11) 


4.  THE  MEASURE  OF  VACANCY 

Define  the  stochastic  process  (I(t),  0<t)  where 


I(t)  - 


,1f  Q(t  mod  2»)>1 
,if  Q(t  mod  2*)«0. 


(4.1) 
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Notice  that  the  sample  functions  of  this  process  are  step  functions,  since  with 
probability  one  there  are  only  finitely  many  random  arcs.  The  measure  of  vacancy 
of  a  specified  arc  from  Ps  to  Ps+t»  t>0,  is  a  r.v.  defined  as 


s+t 

Y(s,t)  *  /  I(t)dT  ,o<s,t<2*.  (4*2) 


This  measure  of  vacancy  is  the  sum  of  lengths  of  all  the  uncovered  arcs  between 
Ps  and  Ps+t*  distribution  of  Y(s,t)  clearly  does  not  depend  on  s  and  is  con¬ 
centrated  on  the  interval  [0,t].  The  corresponding  c.d.f.  is  continuous  on  (0,t) 
and  has  jump  points  at  0  and  t.  Furthermore, 


P  {Y(s,t)  -  t}  -  q(t).  (4.3) 


Let  £r(t)  denote  the  r-th  moment  of  Y(s,t).  In  particular. 


t 

*,(t)  »  E(Y(s,t))  •  /  E(I(t ))  dr  -  £t.  (4.4) 


Indeed,  EU(t))  »  P { I (t )  *  1}  ■ 

Furthermore,  for  every  r>2  and  0<Sj<s2<...<sr<2n, 


r 

E  {  II  T  ( S  -f  )  }  *  p(t.,«.«,  tjr_l). 
i*l 


(4.5) 
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Accordingly,  the  r-th  moment  of  Y{s,t),  for  r> 2,  is 


ir{t)  ■  r!  /  •  •  •  /  E  {  H  I(si)}  ds  ...ds  (4.6) 

1-1  J  «* 

S«S1<...<Sr<S+t 


t  T  T-tj 

*  r!5r  /  di  (t-x)e^f j  dt  e^fti)  /  dt, 

o  o  o  2 


r-3 

T  -  l  ti 

T-t1-t2  1-1  r-2 

ex*(t2)  /  •  •  •  /  dtr.2  exp{X4>(tr.2)  +  M»(?  -  l  t-j)}. 

o  o  i«l 


Let  {»r(t);r>l}  be  a  sequence  of  convolutions  defined  recursively  in  the  follow¬ 
ing  manner: 


,if  t<0 
,if  t>0 


(4.7) 


and 


*r(t) 


t 

J 


o 


<'i(T)<'r-l(t-‘r)dT,  r>2. 
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Accordingly,  formula  (4.6)  can  be  expressed  in  the  form 

t 

£r(t)  *  r!cr  /  (t-T )tp  {2*-T)<ir_i(T)dT.  (4.8) 

o 

From  formula  (4.8)  we  obtain  immediately  that 


d  t 

—  5r(t)  *  rUr  /  ^  (2-ir-r  )^r_i (t  )dx  (4.9) 

o 


and 


d2 

£5  Cp(t)  *  rU^j^ir-tHr.^t).  (4.10) 


This  means  that,  for  every  r>2,  C(t)  is  an  increasing  convex  function  of  t  and 
for  r*l  it  is  an  increasing  linear  function  of  t. 

Introduce  the  moment  generating  function  (m.g.f. )  of  Y(s,t) 


C(v,t)  -  l 
r-o 


tr(t),  -«<v<-  , 


(4.11) 


and  let 


♦(v,t)  *  l  vr  ipr(t )  ,  -»<v<®  . 

r«l 


(4.12) 
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tr 

Notice  that  *(t)<  E{X)  and  therefore  l<*J(t)<l/C.  It  follows  that  *r(t)«0(-y), 

r>l  and  therefore  <|>{v,t)  is  convergent  for  all  real  v.  Furthermore,  from  (4.7), 
the  generating  function  ^(v,t)  satisfies  the  renewal  equation. 


^(v,t ) 


^(t)  +  V  /  w  (t)  t(v,t-T)dr. 
0 


Finally,  from  (4.8),  (4.11),  and  (4.12)  we  obtain  the  formula 


(4.13) 


e(v,t)  -  1  +  Ctv  +  Ev  /  (t-T)  ♦  (2*-t)  ♦(Ev,t)  dT  . 

0 


(4.14) 


Generally,  the  distribution  function  of  Y(s,t)  can  be  obtained  from  (4.14).  We 
provide  now  further  development  for  the  case  of  t<*.  In  this  case  *(2*-t)  «  l 

for  all  0<x<t.  Hence,  from  formulae  (4.7)  and  (4.8)  one  obtains  the  recursive 
formula 


•  «r(t)  «  H  f  *  (u)tr.i(t-u)du  ,  r>2  . 

o 


(4.15) 


It  follows  that  the  m.g.f  of  Y(s,t)  satisfies  the  integral  equati 


on 


w 

5(v,t)  «  r(v,t)  +  5v  /  ^(u)c(v,t-u)du. 


(4.16) 


where 


■r(v,t) 


t 


1  +  £v(t  -  /  ^j(u)du) 
o 


(4.17) 
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Let  h(v,t)  be  the  solution  of  (4.16)  for  all  real  v  and  non-negative  t.  Clear¬ 
ly,  for  0<t<*,  h(v,t)  *  £{v,t).  This  solution  can  be  Interpreted  as  the  m.g.f. 
of  the  measure  of  vacancy  of  a  specified  interval  of  length  t  on  the  real  line, 
when  the  coverage  process  is  by  random  intervals  of  length  X,  having  a  c.d.f. 
F(x).  Let  h*{v,u>),  «>£),  be  the  Laplace  transform  of  h(v,t).  From  (4.16)  one 
obtains 


h*(v,u) 


•y*(v,«o) 

1  -  £v¥*(u) 


(4.18) 


where  y*(v,u>)  and  <r*(w)  are  the  Laplace  transforms  of  y(v,t)  and  ♦  (t),  respec¬ 
tively.  Furthermore,  1 


Y*(v,u) 


(1  -  5vT*(u)) 


(4.19) 


Hence, 


h*(v,w) 


1+  - 
u  ur2 


1 

l~£vf*(u) 


(4.20) 


5.  APPLICATION  TO  A  SHADOWING  PROBLEM 

Consider  a  countable  number  of  randomly  distributed  disks  on  the  plane.  For 
any  Borel  set  B  in  the  plane,  the  number  of  disks,  N(B) ,  centered  in  B  is  a  ran¬ 
dom  variable  having  a  Poisson  distribution  with  mean  um{B)  where  v  is  the  aver¬ 
age  number  of  disks  per  unit  area  and  m{B)  is  the  2-dimensional  Lebesgue  measure 
of  B  (area).  We  further  assume  that,  given  N{B)»n,  the  centers  of  these  n  disks 
are  conditionally  independent  and  uniformly  distributed  over  B.  The  diameters 
of  the  disks  are  i.i.d.  random  variables,  Yj,  Y?,...  having  a  common  c.d.f. 
G(y),  o<y<».  Let  C  be  a  circle  in  the  plane  which  does  not  intersect  any  one  of 
the  random  disks  and  whose  center,  0,  is  uncovered.  The  central  projections  on 
C  of  disks  whose  centers  lie  within  C  will  be  called  shadow-arcs.  The  results 
of  the  previous  sections  are  applied  to  determine  properties  of  the  distribu¬ 
tions  of  the  measure  of  vacancy  of  specified  arcs  on  C. 
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Let  (p,0)  denote  the  polar  coordinates  of  the  center  of  a  disk,  with  respect  to 
0.  Thus,  each  random  disk  is  specified  by  a  triplet  of  random  variables 
(p,e,y).  We  consider  only  disks  with  random  parameter  vector  (p,0,y)  in  the 
set 


S  =  <p<l-y/2;  o<e<2ir,  0<y<l} 


N{S)  has  a  Poisson  distribution  with  mean 


1 


E { N {S} }  *  v  J 


2»  i-y/2 

/  /  pdpdedG(y)  *  2»x, 

o  y/2 


where 


1 


*  ■  r  /  (l-y)dG(y)  . 
c  o 


It  is  assumed  that  G(y)  is  absolutely  continuous  with  p.d.f  g (y ) . 
tional  p.d.f.  of  (p,0,y)  within  S  is 


fs(p.e»y) 


wP9(y) 


2ttX 


,(p,6  ,y)cS 


(5.1) 


(5.2) 


(5.3) 


The  condi - 


(5.4) 


Let  X ( p ,0 ,y )  denote  the  length  of  a  shadow-arc  projected  by  a  disk  with  parame¬ 
ters  (p,0,y).  This  length  is  given  by 


X(p,0,y)  c  2  sin-1  (y/2p) 


(5.5) 


Given  N{S)  *  n,  let  Xj,...,Xn  designate  i.i.d.  random  variables  representing  the 
lengths  of  the  n  shadow-arcs.  The  common  distribution  of  these  random  variables 
is  the  c.d.f. ,  F(x),  of  the  random  arcs  discussed  in  the  previous  sections. 
According  to  (5.4)  and  (5.5) 


(5.6) 


F(x)  *  —  J  (  I  pdp)  g (y)dy 

Ao  y 


2sin(x/2) 


v  D(x)  y2 

—  /  (1-y - cot2(x/2))  g (y ) dy , 

2A  o  4 


where 


D(x)  ■  2sin  (x/2)/(l+sin(x/2))  -  l-tan2(— )  (5.7) 

4 


is  the  largest  diameter  of  a  disk  that  can  yield  a  shadow-arc  of  length  x.  That 
is,  if  y  <D(x)  then  y/(2sin  (x/2))  <  l-y/2.  Consider  the  distribution  func¬ 
tions,  related  to  G(x), 


1  D(x) 

Bjt (x )  «  /  y£(l-y)1_£/2g(y)dy,  o<x<t  (5.8) 

c£  o 


where  ' 


C£  * 


/  y£(i-y)1'£/'2g(y)dy 
0 


is  the  normalization  constant.  Notice  that  c0*2A/u.  The  distribution  function 
F(x)  can  then  be  expressed  as 


F(x) 


0  ,x<0 

c  2 
1  2 

Bq(x)  -  —  B2(x)cot2(x/2),o<x<n 


(5.9) 


♦(t)  «  /  [1-B0(x)]dx  +  ~  /  [l-B2(x)]dx  (5.12) 

t  4c0  t 


+  L  (c  -  ir)  +  -L  [cot  (t/2)  B2(t)  +  (t/2)] 
co  4  2co 


Cl 


6.  COMPUTING  THE  MOMENTS  IN  A  SPECIAL  EXAMPLE 

In  the  present  section  we  apply  the  theory  of  the  previous  section  to  develop 
formulae  for  the  determination  of  the  moments  of  the  measure  of  vacancy  in  the  spe¬ 
cial  case  that  G(y)  is  the  uniform  distribution  on  (0,1).  In  this  case  we  obtain 
according  to  (5.7)  and  (5.8)  the  formulae 
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/  [l-B0(x)]dx 
t 

«  /  tan*  £i)  dx  »  \  tan3(~)  -4tan  (^)  +  (*-t)  , 
{  4  3  4  4 

B,(t)  .  l-|u»»(^£)  *  l 


and 


/  [l-B2(x)]dx  ■  itan5(~)  -  ^  tan3(l^) 
^  5  4  3  4 


+  28tan(~)  -7(«-t) 
4 


Furthermore, 


1 

-  cot  (t/2)  B2(t)  -  D2(t)(l-D(t))1/2 


»  tan  -2tan3(^)  ♦  tan5(-^)  . 

4  4  4 


Finally,  since  c0»l/2,  Cj-4/15  and  c2»l/3  one  obtains  by  substituting 
results  in  (5.12)  the  special  formula  for  *(t) 


,  ,  1  .  .  4  ,1-t,.  4  ,,*-t. 

♦(t)  »  -  T  («-t)  +  -  tan  (—)  +  -  tan3(— ) 


(6.1) 


(6.2) 


(6.3) 


(6.4) 


the  above 


(6.5) 


Accordingly,  in  the  case  of  uniformly  distributed  diameters  on  (0,1),  the  expec- 

«  4  4 

ted  length  of  a  shadow  arc  is  E {X}  ■  ♦(0)»  -  --  +  —  ♦  —  *  0.7306,  and  the  first 

O 

moment  of  vacancy  of  an  arc  of  length  t  is  6j(t)  «  t*exp{-0.7306x). 

In  order  to  determine  higher  moments  4r(t),  r>2,  we  have  to  perform  the  con¬ 
volutions  (4.7)  and  (4.8)  recursively,  where 


+j(t) 


0  ,t<0 

exp  {-^a[(~)  -tan  (~)  -  tan3  (—)]},  0<t<* 


1 


,w<t 


(6.6) 


We  present  now  a  polynomial  approximation  to  the  moments  5r(t),  r>2,  for 
0<t<«.  In  this  range  of  t  values  the  moments  will  be  approximated  by  an  analy¬ 
tic  solution  cf  the  recursive  equation. 


6(t)  -  ^(t) 


C(t)  -  r*  /  *(u)6  (t-u),  r>2,  (6.7) 

o  r-i 


where  *(u)  is  a  polynomial  approximating  ^(u),  over  [0,*].  Notice  that  iMu) 

is  an  analytic  function  and  can  be  approximated  by  a  polynomial  of  a  proper  de¬ 
gree.  We  approximate  the  function  (6.6)  by  the  fourth  degree  polynomial 


*(t)  -  2.0486  -  1.69t  +  1. 124t 2  -  .3489t3  +  0.04111**  (6.8) 


The  coefficients  b-j  of  t1  (i«0,...,4)  in  (6.8)  were  determined  by  the  method  of 
least  squares  by  fitting  a  fourth  degree  polynomial  to  33  points  (t-j (t-j ) ) , 

where  t^ix/32,  i«0 . 32.  The  standard  deviation  of  the  residuals  ♦  (t^)  - 

*(ti),  with,28  degrees  of  freedom  is  o». 00911,  with  a  squared-multiple  correla- 
lation  of  R  «.999.  This  is  a  very  high  degree  of  accuracy  in  approximating 

♦  (t)  by  ?{t).  Define  recursively  the  coefficients. 
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.  1-0 
,  i>0 


and  for  each  r>2 


cr,  1  - 


I  bjCr.i.f.j/Cl) 

j»0  J 


»  i»0,...,4(r-l) 


,  i>4(r-l) 


in  which 


bj  . 


coefficients  of  (6.8)  ,j«0,...,4 

0  »j>4. 


Qne  can  prove  then,  by  induction  on  r,  that 


Cr(t)  -  trcr  ^  ,Cr0  t1  ,  r  >  1 

1-0  ri) 


For  small  arcs,  i.e.,  as  t+0,  formula  (6.12)  can  be  simplified,  bj 

A 

^(t)  by  4»(t )*^p^ (0)+t«p^ (0) ,  where  the  derivitave  of  ^(t)  is 


»;it) 


+  (' '  -rrr )  ■ 


0  <  t  <  it. 


cos  (T) 


(6.9) 


(6.10) 


(6.11) 


(6.12) 

approximating 


(6.13) 
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Thus,  as  t+0  we  obtain  the  approximation 


er(t) 


trcr 


r-1 

I 


j-0 


r-l-j  .  j  j 


J  /  s 

(<) 


(6.14) 


where  and  In  Table  1  we  provide  numerical  values  of  the  nor 

malized  moments  Sr(t)/tr,  for  *“*2 . 10  and  the  limiting  value 

lim  Cr(t)/tr  •  q(t)  ,  (6.15) 


for  values  of  t-j»iir/M,  i«4,8 . 64  (M«64).  The  values  of  the  normalized  mo¬ 

ments  for  the  case  of  1*4  were  computed  according  to  (6.14),  with  * (0)*2. 076326 

and  (0)*-<;1(0).  The  normalized  moment  of  order  1  is  6*. 48162  for  all  t.  In 

Table  2  we  present  the  corresponding  standard-deviations  o(ti),  measures  of 
skewnesss,  Yj(ti)  and  kqrtosis  Y2(ti),  where 

TjUO  *  w3(tl)/®3(ti) 

1  (6.16) 
Y2(ti)  ■  w4(t-j ) / o4*  (t -j )  , 


and  M3(ti),  u4(ti)  are  the  third  and  fourth  central  moments.  According  to  Table 
2,  the  distributions  of  the  measure  of  vacancy  are  for  small  arcs  (as  t-*-0)  nega¬ 
tively  skewed  and  sharply  increasing  near  the  right  limit  of  the  interval.  On 
the  other  hand,  as  t  increases  to  *  the  distributions  become  more  symmetric  and 
can  be  approximated  within  (0, t )  by  Pearson's  Type  I  distributions  (see  Johnson 
and  Kotz  [3]). 


Table  1.  The  Normalized  Moments  cr(t)/tr  Determined 
According  to  (6.12)  and  (6.14), 
for  ti*in/64  (i»4,8,... ,64) 


2 

3 

4 

5 

6 

7 

8 

9 

10 

- 

4 

0.4816 

0.4343 

0.4267 

0.4212 

0.4173 

0.4145 

0.4123 

0.4105 

0.4091 

0.3958 

8 

0.4301 

0.4024 

0.3842 

0.3707 

0.3599 

0.3508 

0.3428 

0.3356 

0.3291 

0.3252 

12 

0.4117 

0.3755 

0.3525 

0.3362 

0.3236 

0.3133 

0.3046 

0.2970 

0.2901 

0.2672 

16 

0.3957 

0.3520 

0.3250 

0.3062 

0.2921 

0.2809 

0.2715 

0.2636 

0.2565 

0.2196 

20 

0.3817 

0.3315 

0.3010 

0.2801 

0.2647 

0.2527 

0.2429 

0.2346 

0.2275 

0.1804 

24 

0.3695 

0.3136 

0.2800 

0.2574 

0.2409 

0.2282 

0.2180 

0.2096 

0.2023 

0.1483 

28 

0.3589 

0.2980 

0.2617 

0.2375 

0.2201 

0.2069 

0.1964 

0.1878 

0.1805 

0.1218 

32 

0.3497 

0.2843 

0.2457 

0.2201 

0.2020 

0.1883 

0.1775 

0.1688 

0.1615 

0.1001 

36 

0.3415 

0.2723' 

0.2315 

0.2049 

0.1860 

0.1720 

0.1611 

0.1522 

0.1449 

0.0823 

40 

0.3343 

0.2616 

0.2191 

0.1914 

0.1720 

0.1577 

0.1466 

0.1378 

0.1305 

0.0676 

44 

0.3280 

0.2522 

0.2081 

0.1795 

0.1597 

0.1451 

0.1339 

0.1251 

0.1178 

0.0556 

48 

0.3223 

0.2438 

0.1983 

0.1690 

0.1488 

0.1341 

0.1228 

0.1139 

0. 1068 

0.0456 

52 

0.3172 

0.2363 

0.1895 

0.1596 

0.1391 

0.1242 

0.1129 

0.1041 

0.0970 

0.0375 

56 

0.3125 

0.2295 

0.1817 

0.1512 

0.1304 

0.1154 

0.1042 

0.0954 

0.0884 

0.0308 

60 

0.3083 

0.2233 

0.1746 

0.1437 

0.1227 

0.1077 

0.0964 

0.0877 

0.0808 

0.0253 

64 

0.3045 

0.2178 

0.1682 

0.1369 

0.1158 

0.1007 

0.0895 

0.0809 

0.0741 

0.0208 
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Table  2.  Standard  Deviations  and  Measures  of  Skewness  and  Kurtosis 


i 

°i 

*1.1 

*2.1 

4 

0.09811 

-0.30546 

1.58630 

8 

0.17480 

0.05021 

1.17473 

12 

0.24975 

0.05300 

1.26210 

16 

0.31780 

0.05611 

1.34884 

20 

0.37993 

0.05902 

1.43394 

24 

0.43699 

0.06130 

1.51655 

28 

0.48976 

0.06262 

1.59601 

32 

0.53389 

0.06290 

1.67165 

36 

0.58490 

0.06221 

1.74299 

40 

0.62823 

0.06089 

1.80958 

44 

0.66921 

0.05942 

1.87096 

48 

0.70808 

0.05835 

1.92677 

52 

0.74505 

0.05814 

1.97704 

56 

0.78027 

0.05897 

2.02223 

60 

0.81398 

0.06057 

2.06333 

64 

0.84627 

0.06179 

2.10328 

20 
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